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SYSTEM AND ME THOD FOR DET ERMINING 
PROPERTIES OF A TUBULAR CAVITY 

FIELD OF THE INVENTION 



The present invention relates to determining properties of a tubular cavity 
including the location of a device inserted into the cavity. 

BACKGROUND OF THE itsVENTION 

In many fields there is a need to determine information about the cavity of 
tubular structures, such as blood vessels or pipelines. This Information can 
include the thickness of the walls, the maximum and minimum internal diameters, 
and the location of a devbe inserted into the cavity relative to the cavity. 

In anatomy the cavities of such tubular organs (for example veins and 
arteries) are called lumens. A typical prior art solution for analyzing lumen 
geometry and navigating within them has been to use ultasound. For example, 
in an ultrasound coronary investigation a full circle scan (2D) is used in order to 
estimate the lumen size of a blood vessel. To achieve the estimate requires a 
large amount of computational power and involves sophisticated image 
processing schemes. 

The present invention proposes a simple method for determining lumen 
properties by analyzing signals from a device placed inside the lumen. 

In summary, it is a main object of the present invention to provide a 
method for analyzing signals from a device inside a tubular cavity to determine 
properties of the cavity including the location of the device. 
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The calculating power required for the present invention is significantly 
less that that required for conventional image processing methods.generally two 
orders of magnitude less. 

Still further objects and advantages will become apparent from a 
5 consideration of the ensuing description and drawings. 

BRIEF DESCRIPTION OF THE INVENTION 

There is thus provided, in accordance with a preferred embodnent of the 
present invention, A computer based method for determining the distance of a 
transceiver located within a lumen from the center of the lumen and for 

1 0 detemnining the radius of the lumen, the lumen crosssection being 

substantially circular at the transceiver location, the method applied on data 
received from a transceiver placed at a position within the lumen that is 
distance (r) from the center and distance (a) firom the lumen wall, transmitting 
a signal of known velocity (v) that can be comelsled with the time of flight and 

1 5 receiving a first signal and a second signal that are reflections of the 

transmitted signal, timing the time differences between the transmission of the 
transmitted signal and reception of the first (t1) and second (t2) refection 
signals, the method comprising: 

a. Calculating the distance of the transceiver fi-om the center of the 
20 lumen = (t1-t2)v/4; 

b. Calculating the diameter of the lumen = (t1 + t2)v/2; 
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BRIEF DESCRIPTION OF THE FIGURES 

The invention is described fierein, by of example only, with reference 
to the accompanying Figures, in which lil<e components are designated by lil<e 
reference numerals. 

FIG. 1 illustrates a model of the geometry of a lumen with a transceiver inside 
it. 

FIG. 2 illustrates a model of the geometiy of a lumen with a transceiver inside 
it taking into account angular reflections along the reflection direction 



DETAILED DESCRIPTION OF THE INVENTION 

10 The present invention proposes a simple method for detemiining lumen 

properties by analyzing signals fran a device placed inside the lumen. The 
device is assumed to be a transducer and the signals are assumed to be 
ultrasound. However the signals can be any signal that can be correlated with its 
time of flight. 

16 Ultrasound measurements consist of an emission of a sound at very high 

frequency, listening to the echo, and analyzing the echo structure. Blood vessel 
wall reflects ultrasound, thus enabling good distance measurement. Various 
components of the blood system reflect ultrasound waves as well, but with 
different efficiencies, i.e. intensities, noise, etc. The sum of those differences 

20 provide information about the blood vessel and the lumen. 

Ultrasound measurements comprise a simple procedure: 

- Emission of a short pulse of acoustic energy 

- Exact registration of the emission time 

- "Listening" to echoes (multiple reflections) 

25 - Electronic and mathematical signal manipulations 
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- Echo classification 

- Distance (anival time) estimates 

- Display in temis of image and/or numbers. 

Ultrasound Is basically a pressure wave or wave packet and behaves 
5 according to simple laws of optics. We do not refer here to nonlinear effects like 
cavitations, frequency dispersion, etc. 

Consequently, we can define our physical system as comprising a pulse 
source (ultrasound emitter), a nondisperslve propagation medium, a cylindrical 
reflector and a receiver (ultrasound receiver). The reflector has two types of 
10 influence on the impinging waves: 

- Reflection; like any other surface it reflects according to the geometry 
and impinging angles. 

- Spatial dispersion: some energy will be reflected around the reflection 
angle, thus a Lambertian element will be present In the reflection 

15 pattern. (A Lambertian reflector is characterized by reflection intensity 

different firom zero at angles oflier than reflection angle. The reflection 
pattern resembles a lobe shape oriented along the reflection direction. 
The more is the reflector Lambertian in nature, the more wide the pear 
shape will be. A perfect Lambertlam reflector reflects light with equal 

20 Intensity at all angles disregarding the direction of the impinging light.) 

The system has a cylindrical symmetry, hence a plane model will suffice. 

The mathematical equations are developed according to the following 
assumptions: 

- Ultrasound source is a point and has a constant angular radiation 
25 pattern. 

- Ultrasound receiver Is a point and has a constant angular sensitivity 
pattern. 
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- Reflector is a circle with a perfectly reflecting mirroFlike sur^ce. 

- Tliere is no intensity attenuation 

The mathematical model is indeed very simple as only two algebraic 
equations are used. The geometry of the model is shown in FIG.1. Lumen (circle) 
5 10 has a radius (R). Transceiver 12 (comprising a transmitter and a receiver, 
both having no effective spatial extension) is located at a distance (a) from the 
lumen perimeter and at a distance (r ) fl-om the center. 

Within the framework of the model, waves are reflected firom a specular 
surface according to the laws of geometric optbs. The rBflecting surface has a 

10 circular shape. Waves emerging from the transmitter will bounce fl-om the surface 
at equal angles with respect to the normal at the surface at the same point. Since 
normal lines at circle are passing through the center they are the same as 
diameters. Hence only reflections along the diameter passing through the 
transducer will reach the receiver. Hence, a receiver can detect only waves 

15 coming from the peripheral points located along the line connecting the circle 
center and the detector. This line intersects the circle at two points. Only waves 
from these points will be detected by the receiver. 

One wave is reflected by the near side and is received at a time given by 

',=^ (1) 

20 where v is the sound velocity in the medium. 

The second wave is reflected by the far side and is received at a time is 
given by 



(2) 



Solving the two equations one has the distance between the transducer 
25 (emitter and receiver) a and the reflector radius R as: 
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4 



4 



(4) 



10 



15 



The two numbers given by equations (3) and (4) provide all the needed 
information: 

r - the location of the transducer within the circle. 
R — The circle radius 

Thus, fi-om a very simple measurement one can provide the "navigation" 
map in terms of path size and location on the path. 

An assumption that can be discarded is the specular minor nature of the 
reflector. A more realistic assumption is a surface that has some angular 
reflections along the reflection direction, I.e. it is Lambertiai to some extension. 
This can be taken into account using the model In FIG. 2. The equations (1) and 
(2) must be modified to: 

rcosQ+acosa=R (5) 

rsinG =asina (6) 

After some simple calculations one finds: 

a = ^^R^+r^-2^r^R^cos9 (7) 

Equation (7) reflects the simple fact that for a variation of the angle 9 
between 0 and 180p, a is monotonously increasing. The significance is crucial: 
The arrival time of every reflection is correlated to the reflection geometry. From 
the point of view of the mirror model this means that all the reflections from other 
points not on the center-to-transducer line will reach the transducer at times 
between /, and in an ordered manner. 
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One should mention that a secondary reflection traveling twice the 
diameter can be detected. Time of flight is precisely fourtimes the result in eq. 
(4). The numerical result is independent of the transducer position inside flie 
lumen. 

5 If one assumes a nanrow Lambertian distribution around the reflection 

angles, then only the small angles reflections will significantly contribi* to the 
line widths. 

Another factor to be taken into account is the optical distortion caused by 
large field of view in the case where the receiver and transmitter are not the 
1 0 same mathematical point). This effect will contribute to some line widening. 

Calculated ultrasound signal shows clearly that the main contribution is 
connected to the line width. 

It should be clear that the description of the embodiments and attached 
Figures set forth in this specification serves only for a better understanding of tie 
15 invention, without limiting ite scope. 

It should also be clear that a person skilled in the art, after reading the 
present specification could make adjustments or amendments to the attached 
Figures and above described embodiments that would still be covaed by the 
scope of the present invention. 

20 



